Decoherence Dynamics of Complex Photon States in a Superconducting Circuit 
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Quantum states inevitably decay with time into a probabiUstic mixture of classical states, due 
to their interaction with the environment and measurement instrumentation. We present the first 
measurement of the decoherence dynamics of complex photon states in a condensed-matter system. 
By controUably preparing a number of distinct, quantum-superposed photon states in a supercon- 
ducting microwave resonator, we show that the subsequent decay dynamics can be quantitatively 
described by taking into account only two distinct decay channels, energy relaxation and dephasing. 
Our ability to prepare specific initial quantum states allows us to measure the evolution of speciflc 
elements in the quantum density matrix, in a very detailed manner that can be compared with 
theory. 

PACS numbers: 03.65.Yz, 85.25.Cp, 03.67.Lx 
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Quantum coherence provides the fundamental distinc- 
tion between quantum and classical states, and is at the 
root of much of the non-intuitive behavior of quantum 
systems; quantum coherence also provides the central im- 
petus for trying to build a quantum computer. Coher- 
ence is however extraordinarily delicate, decaying as a 
quantum system interacts with its environment, or van- 
ishing as the result of even a simple measurement. A 
quantitative understanding of the decoherence process is 
critically important for predicting the behavior of quan- 
tum systems. The harmonic oscillator is a particularly 
compelling system for such studies, as an oscillator in- 
cludes an infinite number of energy levels, from which 
arbitrarily complex states may be created, yet theories 
for decoherence in this system are especially simple. Har- 
monic oscillators are also generic, as they map into a large 
number of different physical systems: A highly useful 
example is the electromagnetic resonator, which is cen- 
tral to a number of approaches to building a quantum 
computer H, H, 0, [H , serving as a quantum memory 
element or a communication bus. The recent develop- 
ment of techniques for generating and measuring arbi- 
trarily complex quantum states in a harmonic oscillator 
@, 0, & B [iqI . [Tl] presents the intriguing possibility of 
exploring the detailed time evolution of these complex 
states, allowin g te sting of the theories of decay dynamics 



Decoherence of complex harmonic-oscillator states has 
been studied in previous experiments [ll, Q • The ef- 
fect of energy decay was measured using Schrodinger cat 
states by observing loss of interference fringes H, 1^ or 
by directly filming the Wigner function |l4l]. Dephasing 
was measured by injecting noise and then observing the 
loss of interference fringes for specially prepared states 
. Here we describe a series of experiments in which we 
prepare specific, complex, quantum-superposed photon 
states in an electromagnetic resonator, and then monitor 
the decay dynamics of its density matrix. By selectively 



preparing certain types of states, we can highlight the 
decay of various elements in the density matrix, in par- 
ticular off-diagonal ones describing phase coherence, for 
a precise comparison with theory. For the first time, the 
time decay of various off-diagonal elements are explicitly 
measured and compared to theory, which show a unique 
and simple signature. We also demonstrate that energy 
decay and dephasing can be measured simultaneously: 
Surprisingly, dephasing rates can be extracted even when 
it is not dominant, indicating the precision of our mea- 
surement. In characterizing a superconducting resonator 
circuit, our experiment measures an unexpected dephas- 
ing rate that is 30 times slower than energy decay. To 
perform these experiments, we have measured the decay 
of more than 13 unique off-diagonal elements, for the su- 
perposed photon Fock states |0) + |n), n = 1 to 8, for the 
states \m) + |3), m = 1,2, and for the Schrodinger cat 
states. 

Our experiment is performed using a half-wavelength 
coplanar waveguide resonator made from superconduct- 
ing aluminum 0, [ll|, • The cavity has a resonance 
frequency of 6.971 GHz, and is weakly coupled through a 
capacitor to an external microwave source. Non-classical 
photon states are generated and measured via a sec- 
ond weak capacitive coupling to a superconducting phase 
qubit, whose state can be manipulated using a second mi- 
crowave source, and measured quickly with single-shot fi- 
delity near unity (~ 90%) [19^. The strong non-linearity 
of the qubit allows complete control of the resonator pho- 
ton state The qubit has an energy relaxation time 
Tiq ~ 300 ns and a phase coherence time T2q ~ 120 ns, 
whereas the resonator has an energy relaxation time 
Ti ~ 2.4 /is, similar to previous devices I2I. [2q|. 

A single measurement of the qubit state gives an out- 
come of 1^) or |e); by repeating the state preparation and 
measurement, typically hundreds of times, we determine 
the probability Pg for the qubit excited state. Tomogra- 
phy of the resonator state involves measuring this prob- 
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FIG. 1: Decay of the |0) + \n) states, (a) Measured density 
matrices at five delay times for the decay of |0) + |3). The 
amplitude and phase of pmn is represented by the length and 
direction of an arrow in the complex plane (for scale, see leg- 
end on right). The element pos (= Pso) is the only off-diagonal 
term that is non-zero during evolution, as expected, (b) Re- 
construction of the Wigner distribution W(a) from data in a. 
The distribution decays to the ground state |0) for long times. 

(c) Amplitudes of the density matrix elements versus time for 
the decay of |0) + |3). Solid lines are fits to the data for pss 
and y9o3, from which decay times T33 and T03 are obtained. 

(d) Plot of Tnn (blue circles), Ton (red squares), and Ton,^ 
(black hexagons) versus n for n = 1 to 8, obtained for a range 
of states |0) + |n). Dephasing times Ton,^ are calculated from 
Eq. [21 Blue line (slope = -1) and black line (slope = -2) are 
weighted least-square fits to obtain the resonator Ti and T^, 
respectively. 



ability for a range of resonator preparation parameters, 
followed by a photon swap with the qubit, as described 
below. 

Non-classical resonator states are generated with qubit 
pulse sequences similar to those in previous experiments 
which are based on the Law/Eberly protocol [l'^. 
The resonator state is subsequently measured after a vari- 
able delay time < t < 15 /is by first displacing the res- 
onator with a coherent microwave pulse, which is charac- 
terized by a complex amplitude a, where |ap corresponds 
to the average photon number. The qubit is then brought 
into resonance with the resonator for a variable interac- 
tion time r, allowing photon exchange between the two, 
after which the qubit state is measured. Because the 
resonator-qubit swapping frequency depends on photon 
number n = 0,l,2...,a measurement of the excited qubit 
state Pg versus r is used to determine the photon occu- 
pation probabilities Pn(<^) in the resonator. In a prior 
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FIG. 2: Snapshots of the time evolution of |2) + |3). (a) 
Measured and simulated density matrices at five delay times. 
The simulation starts with the measured density matrix at 
minimum delay and uses the decay parameters Ti = 2.4 /xs 
and = 70 /xs. Note that only the off-diagonal elements p23 
pi2 poi are non-zero, (b) Wigner distributions VK(a) as 
a function of time, reconstructed from data in (a). 



publication fH] we measured the full Wigner distribu- 
tion W{a) from the parity ^n{—^)^Pn{<^) at thousands 
of points in the complex phase space a. Here we instead 
use a much more efficient sampling for a, of only 60 points 
in total, arranged as two concentric circles with radii \a\ 
= 1.10 and 1.45, from which we calculate the density ma- 
trix and reconstruct the Wigner distribution [21]. The 
resulting readout fidelity is close to that obtained with 
full Wigner tomography, with the increased uncertainty 
(~2%) mostly coming from the smaller sample size. 

The decay dynamics of a harmonic oscillator quantum 
state are described by a Markovian master equation that 
assumes uncorrelated energy relaxation and dephasing 
processes [l8|. The density matrix element pmn in the 
Fock state basis obeys 



dpn 



dt 



m-\-n (m — n)^ 
, 2Ti ^ 
V'(m + l)(n + 1) 



Pm+l,n+l 7 



(1) 



where Ti and are the resonator energy relaxation and 
dephasing times, respectively. For m = n, Eq. [1] reduces 
to the master equation for photon number (Fock) state 
decay [l2j, with the n-photon Fock state lifetime given 
by Tnn = Ti/n. For the off-diagonal elements m ^ 
dephasing causes p^n to decay at an additional rate pro- 
portional to the square of the distance from the diagonal. 
In both cases, decay proceeds along the diagonal pmn 

Pm-l,n-l Pm-2,n-2 ^ ' ' ' (sCC Fig. [S]). 

The decoherence for the Fock state superpositions |0) + 



3 



|n) is particularly simple, as only one independent oflF- 
diagonal element is non-zero. This element, pon (= Pno)^ 
decays at the rate 
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where Ton,0 = T^j^/n'^. By monitoring pnn and pon foi" this 
type of state, we can directly determine the dephasing 
time T(f) and validate the prediction of Eq. O 

In Fig. [T]we show the evolution of the state |0) + |3) at 
five delay times; a full movie of 30 time steps is available 
online [22]. Initially, the only non-zero terms in the den- 
sity matrix are poo, and poa (= Pso)^ as expected. As 
time evolves, the diagonal elements decay to poo- The off- 
diagonal element po3 also decays with time, but without 
other off-diagonal states becoming occupied, as expected. 
In Fig. [Tfc) we display the amplitudes of the P33 and pos 
elements versus time, together with P22, Pu and poo- For 
clarity, statistical errors (~2%) are not shown. From the 
solid line fits we obtain the decay times T33 and T03, and 
using Eq. [2] we extract the dephasing time To3,0. 

We have measured the evolution for the states |0) + |n), 
n = 1 to 8; in Fig. [ijd), we show the measured decay 
times Tnn, Tbn, and Ton,(/). We find Ti = 2.4 ps, with 
Tnn scaling as 1/n, consistent with the n-photon Fock 
state lifetime, as observed previously [HI, \m. We also 
find that Ton,(f) has scaling consistent with confirm- 
ing the dephasing mechanism. The resonator dephasing 
time Tff) is very long, approximately 70 ps, and is roughly 
consistent with off-resonant coupling to qubit dephasing. 
Although Ti and energy decay dominates deco- 

herence, the precision of our measurement allows us to 
indent ify a much smaller dephasing rate. 

Decoherence theory can be further tested by creat- 
ing states with non-zero elements closer to the diagonal, 
which then decay to other elements in the density matrix. 
In Fig. [2] we show the evolution for the state |2) + |3), 
which initially has one non-zero off-diagonal element at 
P23 (= ^32)- the state evolves, this element decreases 
in amplitude, but generates non-zero amplitudes in the 
sequence of elements P23 P12 Poi • 

We plot the magnitude of these elements versus time 
in Fig. [3l along with sequence data taken from the states 
|1) + 1 3) and |0) + |3). The evolution is qualitatively 
similar to that found for the diagonal elements for Fock 
states [HI. We see good agreement between experiment 
and theory, indicating that the off-diagonal evolution is 
as predicted by the master equation. This is the first ob- 
servation of the sequential movement of the off-diagonal 
coherence element in the density matrix due to energy 
relaxation. 

We have also monitored the evolution of states involv- 
ing more than a single pair of off-diagonal density ma- 
trix elements. Fig. [D^a) shows the evolution of the state 
|0) + z |2) + |4), which initially has three unique non-zero 
off-diagonal elements. Fig. St^b) shows the decay of the 
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FIG. 3: Evolution of the off-diagonal elements for \m) + |3) 
with (a) m = 0, (b) m = 1, and (c) m = 2. Right panels 
illustrate the evolution of the off-diagonal elements. Left pan- 
els display amplitude of matrix elements versus time. Solid 
lines are simulations based on Eq. [1] using Ti = 2.4 ps and 
= 70ps. 



coherent state \a = a/5), created by driving the resonator 
with a classical pulse to an average photon number of 
\af = 5 [l3, HH. During state evolution, the coherent 
peak is observed to move toward the origin at a rate con- 
sistent with Ti, while the Wigner distribution remains 
that of a coherent state, as expected from theory with 

Figure [4] also shows the decay of the odd (c) and even 
(d) Schrodinger cat states [14]. These states are cre- 
ated from a superposition of coherent states \a = iV^) ± 
I a = — za/2), where the even (odd) cat is for the plus (mi- 
nus) sign, and consists of a superposition of even (odd) 
number photon states in the Fock-state basis. In the 
experiment, the photon numbers are truncated above 
n = 6. As seen in the figure, the cat states feature a 
Wigner distribution with two peaks, corresponding to the 
coherent states displaced from the origin, but with quan- 
tum interference fringes between them. The even and 
odd cats have opposite fringe patterns, as expected. The 
interference fringes disappear dit ^ 1 ps, indicating the 
loss of phase coherence mostly due to energy relaxation. 
The two coherent peaks begin to merge at ~ 3 ps, indi- 
cating a complete energy decay of the system towards a 
vacuum state. 

Movies of the decay dynamics for all the states men- 
tioned here are available online [22[. In all these states, 
movies constructed from the experimental data agree 
well with simulations based on the measured initial den- 
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FIG. 4: Evolution of (a) the state |0) + i |2) + |4), (b) the 
coherent state \a = a/5), (c) the odd Schrodinger cat state 
I a — iV^) — |a = and (d) the even Schrodinger cat 

state \a = iV^) + |a = — iy^)- Movies of the evolution for 
all states can be viewed online 123]. 



sity matrix and the decay parameters Ti = 2.4 /is and 
T^ = 70/is. 

In conclusion, we have measured the evolution of a 
number of non-classical photon states, prepared and mea- 
sured in a microwave electromagnetic resonator using 
a superconducting phase qubit. The deterministically- 
generated states have non-zero off-diagonal elements in 
their density matrices, and we find that their time evolu- 
tion is in excellent agreement with theoretical predictions 
based on the Markovian master equation, both in the se- 
quence of elements and the individual decay rates. While 
decay from energy loss is in good quantitative agreement 
with theory, we can only conclude that the decay due to 
dephasing is consistent with theory, due to the compara- 
tively small effect of dephasing in this system (T^ ^ Ti). 



We also find that movies of decoherence for all the states 
we prepared and monitored are in excellent agreement 
with simulations. 
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